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We consider the Einstein constraints on asymptotically euclidean manifolds M of dimension n > 3 
with sources of both scaled and unsealed types. We extend to asymptotically euclidean manifolds 
the constructive method of proof of existence. We also treat discontinuous scaled sources. In the 
last section we obtain new results in the case of non-constant mean curvature. 
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Q>^ ■ I. INTRODUCTION 

a^: 

I The geometric initial data for the n + 1 dimensional Einstein equations are a properly riemannian metric g and 

'"^ • a symmetric 2-tensor K on an n-dimensional smooth manifold M. These data must satisfy the constraints, which 
^^ ' are the Gauss-Codazzi equations linking the metric g induced on M by the spacetime metric g with the extrinsic 
curvature iiT of M as a submanifold imbedded in the spacetime {V, g) and the value on M of the Ricci tensor of g. 
As equations on M , these constraints read 

^ . R{g) - K.K + (trK)^ = 2p hamiltonian constraint (1) 

'*' , V.A' — VtrK= 7 momentum constraint (2) 

a\. 

^rt Rid) is the scalar curvature and a dot denotes a product defined by the metric g. The quantity p is a scalar and j 

f.—^ . a vector on M determined by the stress energy tensor of the sources. In coordinates adapted to the problem, where 
Q>^ ' the equation of M in F is a;° = 0, one has 

^: 3^^m\ p=N'T"° (3) 

I ', with p > if the sources satisfy the weak energy condition and if p > g(j, i)^^^ the sources satisfy the dominant 
^ • energy condition. The space scalar N is the spacetime lapse function. 

rr' A classical method of solving the constraints, initiated by Lichnerowicz when n = 3, is the conformal method (cf. m 

and references therein anterior to 1980, 0). In these papers solutions were obtained under the condition that the 
initial submanifold will have constant mean extrinsic curvature, i.e., trK = constant. Recently the results have been 
. ■ . extended to the non-constant mean curvature case with some hypotheses on the smallness of its variations. The case 
C^ ' of a compact manifold M is treated in 1^ and 0, the first by using the Leray-Schauder theory, the second through a 
constructive method. Results for asymptotically euclidean M are given in M , using again the Leray-Schauder theory. 
All the quoted papers treat the case of scaled and continuous sources on a three-dimensional manifold M . 

We will in this article consider the case where the manifold M has an arbitrary dimension n > 3 and the sources are 
the sum of scaled and unsealed ones. We will extend to asymptotically euclidean manifolds the constructive method. 
We will extend the existence proof to discontinuous scaled sources. 

In the last section we obtain results in the non-constant trK case. In the asymptotically euclidean case, non-constant 
trK denotes non-maximal submanifolds. A simple smallness assumption on the variations of trK is sufficient to insure 
existence of solutions for metrics in the positive Yamabe-Brill-Cantor class when there are no unsealed sources. In 
the other cases the study is more delicate, as pointed out by O'Murchadha, and we obtain some results, in particular 
for unsealed sources. 

We do not claim to have constructed solutions with scaled sources in the negative Yamabe class on non-maximal 
manifolds. The problem of the existence of solutions with large variations of trK remains also open. 

We will use the conformal thin sandwich formulation developed recently by one of us [|6| to express the momentum 
constraint. It gives a better understanding of the splitting between given and unknown initial data. 
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II. CONFORMAL METHOD IN ITS THIN SANDWICH FORMULATION 

One turns the haniiltonian constraint into an elliptic equation for a scalar function ip by considering the metric g 
as given up to a confornial factor. A convenient choice is to set, when n > 2, 

2 

g = l"/^ , «-e. mj = <p^^7*j , withp = ——^ . (4) 

Then the following identity holds: 

Rij9) ^ ^-("+2)/("-) (^^i?(^) - ^^A,^) . (5) 

The haniiltonian constraint becomes a semi-linear elliptic equation for Lp with a non-linearity of a fairly simple type 
when 7 and K are known — namely 

A^(^ - fc„i?(7)(^ -f K[K.K - r2 + 2p)(^("+2)/("-2) ^ q (g) 

with 

r = trK, fc„^ "~^ . (7) 

4(n — 1) 

We now explain the conformal form of the momentum constraint as recently deduced by one of us [pj from thin 
sandwich considerations. It can be construed to include previous methods as special cases, but no tensor splitting is 
needed. The initial metric g being known up to a conformal factor, it is natural to consider that the time derivative 
of this metric (the other ingredient of the initial data in a thin sandwich formulation) is known only for its conformal 
equivalence class. We had above 

m, = ^'/("-'^7., • (8) 

If gij and 7^ depend on t, their time derivatives are linked by 

with 

dtgtj 9ij9^''dtghk = Uij (10) 

and an analogous expression for Uy constructed with 7^ . 

We will consider the traceless symmetric two-tensor u^ as given on the manifold (A/, 7) . Recall the identity 

K,, = {2N)-^{-dtg,j + V J, + V,h} , (11) 

where (3 and N will be respectively the shift and the lapse in the imbedding space time. The shift vector /?' is not to 
be weighted; it is not a dynamical variable. The other non-dynamical variable is not the lapse N but a scalar density 
a of weight —1 such that N = adet(5)^/^ (cf. Q). We therefore consider as given in this context a function N with 
the space time lapse N linked to it by the relation: 

N = ^2„/(n-2)^ (^2) 

We denote by V and V the covariant derivatives in the metrics g and 7 respectively. We denote by C the conformal 
Killing operator 

{cpy^ = v'p^ + v^p' - -f^VhP^. (13) 

n 
We have 

(£^)^^'ee^-4/("-2)(£/3)^J', P' = (3' (14) 

and 



n' 
with 

A'^ = [2N)-^{-u'^ + (£/3)*^'} . (16) 

One finds by straightforward calculation that the momentum constraint now reads as an equation on (M, 7) with 
unknown /? (and ip if Dt ^ 0) : 

V,{(2iV-i)(/:/?)*^'} = Vj{(2iV-i)u'^} + !^— i^2„/«-2yz^ 

where N , t, and w are given. 

The hamiltonian constraint now reads 
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(^("+2)/(«-2)^2 ^ _2fc„p(p("+2)/("-2). (18) 

4n 
The sources are decomposed into scaled and unsealed sources by setting: 

J = J+^-2(n+2)/(„-2)^^ J^ - ^ p = C + gy,-2(»+l)/(n~2) (^g) 

2(n — 1) 

More refined decompositions may also occur (See example 2 below). 

The energy density scalar p and the momentum density vector j behave under conformal rescaling of the metric 
according to the source fields which they represent. See references by Isenberg, O'Murchadha, and York and by 
Isenberg and Nester in y . 

Examples. 

1. n = 3, the source is an electromagnetic or Yang-Mills field F. The electric and magnetic fields relative to a 
spacetime observer at rest with respect to the initial manifold M, (i.e., with 4- velocity orthogonal to this manifold), 
are 

Ei ^ N~^F'q = Lp-^N-^Fl = tp-^E' (20) 

H' = Iv^'^F.k - ^^-N^^'=F,fc = ^-''W (21) 

with 77 and fj respectively the volume forms of 7 and g. 

Note that if {E\H^) satisfy the MaxweU constraints V^i?* = and V^if'' = in the metric g, the fields {E\W) 
satisfy these constraints in the metric 7. We consider that it is these last fields which are known on M. 

The energy density is 

p ^ \g,j(&E^ + H'H^) = ^-\ (22) 

with (7, considered as known on M, given by 

q=]p,,{E'E^+WW) (23) 

The momentum density is 

■i ^ J^J.^o ^ ^^Oj^j ^ -E'g'^T^kjiH^ = ^^^^v' (24) 

with u* the quantity considered as given; 

v' = -Y^VkjiE^H' . (25) 



The sources are scaled as defined and the constraints decouple if Dt = 0. Note that if g > {'-fijv'^v^)^ , then 

2. General n, the source is a Klein-Gordon field. The energy density on M of a Klein-Gordon field ip with respect 
to an observer at rest is 

p = i(iV-2|aoV|' +5'^5,V5,V + mi;^) , (26) 

i.e., 

P = lw^-N-^\do^\^ + ^^7^^9,VajV + mi^^} ■ (27) 

If we consider as known on M the initial data i/'Uf and doiplM together with 7 and N, then neither of the terms in 
p is properly scaled as indicated in (19). The term N~'^\do(p\'^ adds in the Hamiltonian constraint to A. A, the term 
771^-0 is unsealed and gives a contribution to c, the middle term gives a new, positive, contribution to the (p term 
which adds to —R^j). The momentum density is 

f = -N~^g''dj^dotP = -(^-2(n+2)/(n-2)yj^^.^^^^ (28) 

We see that the momentum is properly scaled. The constraints decouple if Dt = 0. 

The methods we give below to study the constraints with properly scaled or unsealed sources can be applied to 
more general scalings, such as this example. 

Summary. The given initial data on a manifold M are on the one hand (geometric initial data) a set (7, u, r, N), 
with 7 a properly Riemannian metric, u a traceless symmetric 2-tensor, r and N scalar functions, and on the other 
hand (source data) a set (J, v, c, q), two vectors and two scalars. The initial data to be determined by the constraints 
make a pair (1^, /3) with ip a scalar function and /3 a vector on M. In the conformal thin-sandwich formalism the 
constraints reduce to the equation (|l^) and (^ which read, taking ( p^ into account, 

W,{{2N-^){C(3y^} = h\.,p) (29) 

with 

h'(., ip) = V,{(2iV-l)w*n -I- !llli-(p2«/(n-2) VV + ^2{n+2)/{n-2)J^ ^ ^^ ^ /ggx 

n 
and 

A^^ = f{.,p), (31) 

where 

/(., ^)=rip- a(p(-3«+2)/(n-2) ^ ^^(„+2)/(„-2) _ q^-n/{n-2) ^ (33) 

with r, a, and d defined as functions of the geometric data as in equation (45). 

When T is constant on M and the sources have no unsealed momentum (i.e. J — 0) these constraints decouple in 
the following sense: the momentum constraint (E9) is a linear equation for /3, independent of (p, and the Hamiltonian 
constraint ( pl| ) is a non-linear equation for (p when f3 is known. 

When the constraints are solved the spacetime metric reads on M: 

ds^ = -N^dt^ + g,j [dx' + (3'dt){dx^ + P'dt) , (33) 



with g and N given by the formulas (|g) and ([121). The extrinsic curvature of M is determined by (y_5|) and (|16|), the 
derivative dtfjij on M by dO). The derivatives dtN and dt(3 remain arbitrary. 

We now express in our setting the conformal invariance of the conformal constraints. 

Lemma The constraint equations (O) and (Eq) are conformally invariant in the following sense: If (/3, ip) is 
a solution of the constraints with data {■y,u,T, N; J,v,c,q) then {I3,ip) is a solution of the constraints with data 

(7 = ((^^-l)4/(«-2)^~ ^ ((^^-l)4/("-2)M,r,7V = ((^(p-l)2«/("-2)^. J ^ J^y ^ (^^-l)-2(n+2)/(„-2)^^~ = C, q = 
((^(/J-l)-2(«+l)/("-2)g). 

Proof. If (/3, ip) together with the considered given data is a solution of the conformal constraints, the corresponding 
Einstein initial data set {g, K) is a solution of the Einstein constraints with sources j, p given by (p^. The Einstein 



initial data set and sources constructed with the ~ quantities are identical with {g, K) and (j, p). Since the Einstein 
constraints are satisfied, the conformal constraints written with the r^ quantities are also satisfied. 

Remark. In the case n = 2, equations analogous to the ones obtained here for the conformal factor ip and the vector 
(3 are obtained by setting (cf. 0): 

.9 = e2^7 , (34) 

and in the thin sandwich point of view, 

N = e^^N (35) 

which gives: 

KV ^ ^^'PA^3 + \gio^ (36) 

However we will not consider n — 2 because it poses special problems in what could correspond to an asymptotically 
euclidean case. 

III. ASYMPTOTICALLY EUCLIDEAN MANIFOLDS AND 
WEIGHTED SOBOLEV SPACES 

The euclidean space E" is the manifold R" endowed with the euclidean metric which reads in canonical coordinates 
^{dx^Y' . A C°°, n-dimensional riemannian manifold (M, e) is called euclidean at infinity if there exists a compact 
subset S of M such that M — S in the disjoint union of a finite number of open sets Ui, and each {Ui, e) is isometric 
to the exterior of a ball in E". Each open set Ui d M \s sometimes called an end of M. If M is diffeomorphic to R", 
it has only one end; and we can then take for e the euclidean metric. 

A riemannian manifold (M, 7) is called asymptotically euclidean if there exists a riemannian manifold (M, e) eu- 
clidean at infinity, and 7 tends to e at infinity in each end. Consider one end U and the canonical coordinates x^ 
in the space E" which contains the exterior of the ball to which U is diffeomorphic. Set r = {X](^*)^}^^^- ^"^ ^^^ 
coordinates x'^ the metric e has components e^ = 6ij. The metric 7 tends to e at infinity if in these coordinates 
•jij — Sij tends to zero. A possible way of making this statement mathematically precise is to use weighted Sobolev 
spaces. (One can also use in these elliptic constraint problems weighted Holder spaces, but they are not well adapted 
to the related evolution problems). 

A weighted Sobolev space W^ ^ , 1 < P < 00, s G N+, 5 G M, of tensors of some given type on the manifold (Af, e) 
euclidean at infinity is the closure of C^ tensors of the given type (C°° tensors with compact support in M) in the 



norm 



\w^. 



i/p 
J2 I |a'"u|P(H-d2)^p(*+™)d/ii , (37) 



a<m<s ^ 



where d, \ \ and dfi denote the covariant derivative, norm and volume element in the metric e, and d is the distance 
in the metric e from a point of M to a fixed point. If (M, e) is a euclidean space one can choose d = r, the euclidean 
distance to the origin. We recall the multiplication and imbedding properties (cf. pJj^ ) 

Tl Ti 

Wl,s, X wis, C Wl, \is<Si,S2, S<Si+S2~- ,5< 61 +62 + -, 

77 77 

W^s C C^ if m<s--, (3<5+-, 
^ p p 

\\u\\c^^ J2 sup(|5^u|(l + d2)i(/3+^)). (38) 

0<e<m ^" 

The imbedding of the space W^g into W^, g,,s>s', 5 > 5' is compact if s > s', (5 > 5' . We have on the other hand: 

[i + d^yf^'^ewi^ a p>5 + - , s>o. (39) 

Let (Af , e) be a manifold euclidean at infinity. Then the riemannian manifold (Af , 7) is said to be 'W^ „ asymptotically 
euclidean" if 7— e G W^ . When we speak of "asymptotically euclidean manifolds" without further specification, we 
suppose that 7 — e G W^p with (T>^-|-l,p>— ^. These hypotheses imply that 7 is C^ and 7 — e tends to zero at 
infinity. 



IV. MOMENTUM CONSTRAINT 

In the thin sandwich confornial formulation the monientuni constraint reads 

W,{{2N-^){Cpy^}^h{.,^) (40) 

with 

where N and t are given functions on M and u a given symmetric traceless tensor field. The sources J and v are 
also considered as known. We suppose momentarily that (p is also a known function; in fact, it disappears from the 
equation if Vr = and J = 0. 

The momentum constraint is a linear elliptic system for the unknown f3 on the manifold (ill, 7). (The symbol of 
the principal operator is an isomorphism.) 

Theorem. Let (M, 7) be a W^ asymptotically euclidean manifold with a > - + 2,p > — — . Let u, t G Wf,-^^ ^.-^ 

be given, (1 — A^^^) and {1 — ip) G M^f+2 51 ^ > Oi '/' > 0, and J,v £ W^g^2- The momentum constraint has one and 
only one solution (3 £ Wf^^^ if s > - - 2 and < s < cr - 2, -^ < (5 < n - 2 - 2. 

Proof. The operator on the left hand side is injective on ^^4.2 5 because a solution /3 € W^g, 5 > — -, of the 
equation 

Vj{{2N)-^{CI3Y^} = Q (42) 

is necessarily a conformal Killing field. Indeed if /3 e C5" the equation implies by integration on M that 

AVj {{2N)-\C(3y^] ^i^ = [ {2N)~^C(3Xpn^ = 0. (43) 

M J M 

The same is true if /3 G Wf^a.S' with (5' > -^ + f - 2 (respectively 5' >~2iip^ 2). There is such a (5' if /3 s W'f+a.i 
satisfies the homogeneous second order equation (cf. a similar proof for the Laplace operator in the appendix). It 
is known that there are no conformal Killing vector fields tending to zero at infinity on an asymptotically euclidean 
manifold (cf. R]) where a proof requiring only low regularity is cited). 

Because the elliptic operator on (3 is injective, the isomorphism theorem applies to give the existence and uniqueness 
of^. 

V. HAMILTONIAN CONSTRAINT 

In the conformal method the hamiltonian constraint reads as a non-linear elliptic equation for the confornial factor 
(f. We write it 

A^t/J = /(.,¥') 

/(., p)=r^- a<^(-3»+2)/("-2) + ^^(«+2)/(„-2) _ q^-n/{n-2) ^ (44) 

with A given by dlq)), and 

r = knR{-f), a = k^A.A, k^ = (n - 2)/4(n - 1) 

d=b-c, b={n-2)/{4:n)T^. (45) 



By their definitions we have 



a > 0, 6 > 0, c > 0, q>0. (46) 



The functions q and c, scaled and unsealed sources, are considered as given on M. We will suppose that r (hence 
b) is also known on M . The function a is known when the momentum constraint has been solved: this can be done 
independently of (^ if r is constant and the unsealed sources have zero momentum. 



The constructive method of sub and super solutions used by one of us ||] to solve non linear elliptic equations on 
a compact manifold can be extended to asymptotically euclidean manifolds. 

The following theorem is a particular case of the theorem proven in the Appendix B. 

Theorem. Let (M, 7) be a {p,a,p) asymptotically euclidean manifold with a > - + 1, p > — — . Suppose r, a, 
q, d £ W^g_^2^'^ — 1 > s > f ,— -^ < S. Suppose the equation A^(p = f{x,ip) admits a subsolution i^_ > and a 
uniformly bounded supersolution (p+ , functions in C^ such that 

A^V'- > /(•, V-) , A^^+ < /(., ^+) (47) 

and 

lim</3_ < 1 , lim(^+ > 1 

V- < ^+ on ^^ ■ (48) 



00 



Suppose that Dip_,Dip+ e VKf,_^ g,_^_^ , s' > s , S' > --. 
Then the equation admits a solution (p such that: 

v-<^<v+, 1 - "^ e Tyf+2.5 (49) 

if 

--<S<n-2--. (50) 

P P 

Remark. When r = knR{j) we have r € W^_2 5,2 ii(J>- + 2,S>~-. We will use this theorem directly in 
Section XI, with constant sub and super solutions. We will give and use in Sections |V| and ^ intermediate simple 



steps to obtain non-constant sub and supcrsolutions. 

VI. BRILL-CANTOR THEOREM 

The constraints in their conformal formulation are invariant under conformal rescaling (cf. Section ||) . 

In the case of a compact manifold M a convenient first step before studying the solution of the Lichnerowicz 
equation is to use the Yamabe theorem which says that each manifold (M, 7) is conformal to a manifold with constant 
scalar curvature which can be chosen to be 1, —1 or zero. The positive, negative and zero Yamabe classes correspond 
to the signs of these constants and are conformal invariants. There is no known analogous theorem for asymptotically 
euclidean manifolds. (In any case the curvatures could not be non-zero constants.) However an interesting theorem 
has been proved by Brill and Cantor, with the following definition. 

Definition. The asymptotically euclidean manifold (M, 7) is in the positive Yamabe class if for every function / on 
M with / e C^ , / ^ 0, it is true that 

/ {\Df\'+ri^)f}p^>0. (51) 

Jai 

The positive Yamabe class is a conformal invariant due to the identity 

A^/ - r(7)/ EE ^("+2)/(-2)|A^,/' _ r{j')f'} (52) 

7' = ^V(«-2)^ ^ /' = /^-i, (53) 

which gives after integration by parts with / e C^, because fi-yi = iy9^"/"~^/i^, 

{\Dff+r{^)f}p,^ f {\Dff+ri^')f"}p,, (54) 

M JM ^ ' 

We will say, following O'Murchada that the asymptotically Euclidean manifold (M, 7) is in the negative Yamabe 
class if it is not in the positive one [[7|. However, analogy with the case of a compact manifold can be misleading, as 
shown in the following theorem. 

Theorem. ( [Q). The asymptotically euclidean manifold (Af, 7) is conformal to a manifold with zero scalar 
curvature, that is, the equation A^(/? — r(pi)ip = has a solution 93 > 0, if and only if (M, 7) is in the positive Yamabe 
class. 



The physical metric g that solves the constraints together with the symmetric two-tensor K has a non-negative 
scalar curvature R{g) if the sources have positive energy and the initial manifold has constant mean extrinsic curvature, 
(necessarily zero in the asymptotically euclidean case).. Thus, R{g) > 0, with R{g) ^ except in vacuum for an 
instant of time symmetry, i.e. K = 0. Therefore, the physical metric g on an initial maximal submanifold is in the 
positive Yamabe class and all metrics 7 used as substrata to obtain it must be in that class. 

We will prove a more general theorem. We will also make fewer restrictions than Brill-Cantor on the weighted 
spaces. 

Theorem. On a (p, a, p) asymptotically euclidean manifold the equation 

A-y(/3 — aip — V , (55) 

where a, v G Wf ^_,_2, v < 0, has a solution (^ > 0, (^9 - 1 G Wf+2 5; * > ; 5 > -- only if for aU / e C^ , / ^ 0, the 
following inequality holds: 

/ {\DJ\'' + af]^l^>Q. (56) 

J M 

Under the same hypothesis the solution ip exists with (p — 1 G Wj 1 2 s^ ^^"-^ <p > if one supposes moreover s > — , 
5 > 2. — s._iifp^2 (respectively 5 > —1 if p = 2), and that either v < or v = on M, or a = r(7) with, in this 
last case ct > 2. 

The theorem of Brill and Cantor corresponds to the case v = and a ~ r(7). They made the additional hypothesis 
p > n. 

Proof. 
1 . ( "only if " ) Suppose ^ exists and solves the equation satisfying the hypothesis of the theorem. Then we will show 
that for any /^O, f e C^ , 

f {\Df\^ + af}^i^>0. (57) 

Jm 

Indeed, let / e C^ , f ^ 0. The function 6 = f(f^^ has compact support, belongs to W2 gi for any S' and is such that 
DO ^ since 6, having compact support, cannot be a constant without being identically zero. We have by elementary 
calculus: 

\Df\^ ^ \D0\^(f^ + ^D(f.D{0^) + 9^\D(f\^ . (58) 

The following integration by parts holds for the considered functions: 

(pDip.D{e^)fi^ = f -e'^D.{LpD(p)ii^ ; (59) 

IM JM 

therefore, 



and 



ipDtp.D{e'')^i^ = -e^itpA^ip+\Dip\^)n^ (60) 

M Jm 



f |i?/|V7 > / ~o^^A.,ip^l^ . (61) 

Jm Jm 



'M JM 

Hence when (p satisfies the given equation and 6ip — f: 

f {\Df\^ + af}ti^> f -v0Va^7>O (62) 

Jm Jm 

a (f > and V < 0. 

2. ("if ":existence) Setting p = 1 + u the equation reads: 

A-yU — au = V + a . (63) 



The operator A^ — a is injective on W^ g (cf. Appendix A). 

The general theorem on linear elliptic equations on an asymptotically euclidean manifold shows that our equation 
has one solution u E VFj^_2 g, s > 0, — - <S<n — 2— -. The problem is to prove that ip ^ 1 + u is positive. We 
will use the maximum principle, supposing the solution to be C^ ,i.e., s > ^. Since a is not necessarily positive we 
cannot apply directly the maximum principle. One proceeds as in the Brill-Cantor proof. One considers the family 
of equations, which all satisfy the criterion for the existence of a solution ip\ with ip\ — I E VFj, 2 g, 

A^Lp - Xaip = Av , A e [0, 1] . (64) 

The solutions ipx depend continuously on A and we have (po ~ I. If the function ifi = ip takes negative values there 
is one of these functions ip\^ which takes positive or zero values. The points where ip\^ vanishes are minima of this 
function. It is incompatible with the equation satisfied by ip\g if v is negative at that point. Therefore we have ip\ > 
for A G [0, 1] if V < 0. 

To prove that (fx„ > 0, and hence ipx > for A e [0,1], when v = 0, we use, as Brill-Cantor, a theorem of 
Alexandrov: if there is a point xq where ip\„ = 0, it is a minimum of this function, hence Dip\g{xo) = 0. Since the 
function ip\g and the function identical to zero take the same value as well as their first derivatives at xq and satisfy 
the same elliptic equation they must coincide (Alexandrov theorem), a result that contradicts the fact that ip\g tends 
to 1 at infinity. 

If we know only that v < but a ~ r{j) wc first conformally transform the metric 7 to a metric 7' — 7^4"/("^2) 
with zero scalar curvature: this is possible by the previous proof for v = (original Brill-Cantor theorem). The 
equation to solve is equivalent to the following equation for ip' = pip^^: 

Ayip' ^ ^-i^'+^)/("-^)y < . (65) 

whose solution is (/?' > 1 because p' cannot attain a minimum at a point of M and (/?' tends to 1 at infinity. 

VII. SOLUTION OF THE EQUATION 



Theorem. If 6 G W^g^^' s > ^,-^ < 5 < n - 2 - ^,b> 0, the equation 

A^ip - r{"f)p = lj^(n+2)/(n-2) (gg) 

on the (p, (7, p) manifold (M, j), a > - + 2, p > - has a solution p = 1 + u,u E W^_^2 51 * < o'j P > under one or the 
other of the following hypotheses: 

1. On the subset of M where r{j) < there exists a number /i > such that 

sup ^^ < ^l . (67) 

{xeM,r{-f){x)<0} 

2. (Mj-f) is in the positive Yamabe class. 
The solution is unique in both cases. 
Proof. 

1. The manifold (M, 7) and the function f{x,y) = r(7)0 + fe((/))("+^'/"^^^ satisfy the hypothesis (H) spelled out in 
Appendix B. The equation admits the subsolution p^ = 0. A number (p_|_ is a supersolution if 

p+ > 1 and r(7) + bpf^'-''^^^ > on M . (68) 

The second inequality is a consequence of the first if r{'j) > 0. 

The hypothesis made on (M, 7) on the subset r{"f) < insures the existence of the number (p_|_ > ip^ = 0, given by 

^+=max(l,^("-2)/4n)_ (gg) 

The existence of a solution -0, with < '0 ^ 'P+ and 1 — ip E Wj 1 2 g results from the general theorem. Such a solution 
can be obtained constructively. We know that ip ^ Q since it tends to 1 at infinity. 



We show that V' > on M by using the Alexandrof theorem as we did in the proof of the Brill-Cantor theorem: if 
■0 vanishes at a point xq ^ M this point is a minimum of ip, hence Dip(xo) = 0. The functions ip = ip and (^9 = both 
satisfy the elliptic equation 

A^(^ - (r(7) + 67A4"/"-2)<^ = . (70) 

They as well as their gradients take the same values, zero, at the point xq, therefore they coincide. This contradicts 
the fact that ip tends to 1 at infinity, therfore there exists no point .tq where ^'(^^o) = 0. Hence ^ > on M . 

2. If (M, 7) is in the positive Yamabe class we conformally transform it to a manifold (M, 7') such that r{^f') = 0. 
The subset of M where r{'^') < is empty; therefore, (/?+ = 1 can be chosen as a supersolution. The proof that (p > Q 
on M can be made using simply the maximum principle: a solution </? G C^ of the equation 

A^Lp - 5^("+2)/(n-2) = A^(^ - (6^4"/("-2))^ ^ Q (71) 

with 6 > cannot attain a nonpositive minimum on M without being a constant (which is not possible with ip tending 
to 1 at infinity except if 6 = 0, in which case (^ = 1). 

The uniqueness property in case 2 is simply a consequence of 6 > and of the increasing property with (^ > of 
the function (^("+2)/("-2)^ together with the fact that the difference of two solutions tends to zero at infinity. The 
uniqueness in the general case results from the conformal properties. Indeed suppose the equation 

A^p-r{^)p^bpQ ,Q='^ (72) 

n — Z 

has two solution pi and p2. We deduce from the conformal identity 

2 

A-fP ~ r{-f)p = -r{g)p'^ , g = p^^-j , p = (73) 



that 



n " I 'T'i /^i I n ' \ 



{g) = K7^f ) - K7^f ) = -h ■ (74) 



Consider the identity 



^i^l^iVi V2) - r{-ipf)p^ V2 = -r{"ipf){Lp^ V2)'^ • (75) 

It implies, because of the previous equalities, 

uQ-i - 1 
'^Ta,2p(w- 1)-^" -. — (w-l) = 0, u = p^^p2. (76) 

We have 6 > 0, w > 0, ^", _■^^ ' > since u > and Q > 1. We deduce from the fact that u — 1 tends to zero at 
infinity that u — - 1 = on M, i.e., (^1 = p)2- 

Remark. By the above theorem, under the hypothesis made, an asymptotically euclidcan manifold [M, 7) is 
conformal to a metric 7' of given non-positive scalar curvature r(7'), and the solution p of the equation (66) with 
b — -~r{'-f') gives the conformal factor. (This result was known to O'Murchadha.) 

VIII. SOLUTION OF THE EQUATION 

Aj(p~R{^)p + A(p-'' + Q(p-'"^0, A>0,Q>0 

This equation is the conformal expression of the Hamiltonian constraint on a maximal manifold with no unsealed 
sources. The following theorem has been proved independently in the case n = 3 in 1979 by Cantor, and Chaljub- 
Simon and Choquct-Bruhat (in weighted Holder spaces). We give here a new constructive proof; the corresponding 
function f{x, p) satisfies the hypothesis H of Appendix B on any interval [£, 00), d > Q. 

The generalized Brill-Cantor theorem shows that the considered equation can have a solution p > Q only if (M, 7) 
is in the positive Yamabe class, a result in agreement with the fact that the original hamiltonian constraint on an 
initial maximal submanifold {M,g) implies r{g) > 0. 

Theorem. The equation on the (p, ct, p) asymptotically euclidean manifold (A'/, 7), ct>- + 2, p>^— 2— -if 
p ^ 2, and p > —1 ii p = 2, given by 
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P == (3n - 2)/(n - 2) , P' ^ n/{n - 2) , 

Tl ft Tl 

a , q e W^f . , „ (7-2>s> - , n-2- - > S> , (77) 

s,o-|-z p p p 

has a solution ip > 0, tp — 1 E W^,^ g if and only if (A'/, 7) is in the positive Yamabe class. This solution is such that 
(p > 1. It can be obtained constructively and is unique. 
Proof. 

1. ("only if) This part follows from the generalized Brill-Cantor theorem. 

2. ("if) The manifold (M, 7) is conformal to a manifold (M, 7') with zero scalar curvature, 7' = ^4/("~2)^ ^ r(7') = 0. 
Conformal covariance shows that the resolution of the given equation is therefore equivalent to the resolution of an 
equation of the same type but with no linear term, which, suppressing primes, we write as 

/\~,ip = -a^-^ - qip-^' , a > , g > . (78) 

This equation admits a constant subsolution ip- = 1 but no finite constant supersolution. However, it is possible to 
construct a sequence u^ G ^^+2 .5 starting from the subsolution </?_ — 1 by solving the equations with A: > , A: G 

AjU,y — kui, = — a(l + Uy_i)^-^ — q{l + u^_i)^-^ — ku^^i . (79) 

We have Ui, € VFj_^_2 s '^ ^a ^'^'^ ^11 '^ such that a < S + -, hence u^ tends to zero at infinity and we can use the 
maximum principle to see that Ui, > 0. We could choose k > Pa + P'q and deduce as before from the maximum 
principle that the sequence Ui, is pointwise increasing, but we do not obtain an upper bound through the maximum 
principle because we do not have a supersolution. We choose first instead fc = to construct our sequence and write 
the elliptic estimate, using the fact that (1 + u^)^^ < 1 since u^ > , 

\K\\wl,<c{\\a\\^.^^^ + \\q\\^.^^^} . (80) 

The sequence, being uniformly bounded in T4^f ^, admits a subsequence which converges in the W^ g, norm, S' < S, to 
an element u G W^f 5- The rest of the proof is the same as in the general arguments given in Appendix B, except that 
in the present case, the sequence Ui, is not proven to be monotonic, nor identical to the subsequence which converges. 
Hence we cannot conclude that the limit u of the subsequence is a solution of (78). 

To obtain a converging sequence, and consequently a solution, we again use (79), but now with k > Pa + P'q. 
For (79) with such a k, the subsequence limit u serves as a supersolution. Therefore, the increasing sequence Ui^ is 
bounded above by u and it converges to it in W^,^ 5 ■ We have ip > I. A pointwise upper bound for (p can be deduced 
from the Wj, 2 s norm oi u = ip — 1 . 

Remark. The sequence Ui, and the limit u, bounded in W2 g norm in terms of the W^ g,2 norms of a and q, are 
therefore bounded in C° norm in terms of these norms of a and q if p > §. 

3. Uniqueness: the equation with r(7) = has a unique solution such that ip tends to 1 at infinity because of the 
nionotonicity of the right hand side and the maximum principle. The original equation also has a unique solution. 

IX. SOLUTION FOR SCALED SOURCES 

We now prove an existence theorem for the non-linear elliptic equation for (p expressing the Hamiltonian constraint 
on an arbitrary initial manifold, when there are no unsealed sources. 
Theorem (scaled sources). The equation: 

A^ip - r(7)^ = /(., ^) = -a^-P - q^-^' + 6^^ (81) 

with a> 0, q> 0, b> 0; a,q,b e W^ gj^^.^ s>-,--<S<n-2--, has a solution ip = 1 + u, u e WJY2 5' "^ > 
which can be obtained constructively, if either (a.) or (b.) holds: 



(a.) On the subset where r{"f) < 



|r(7)|<6, (82) 
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(b.) (M, 7) is in the positive Yamabe class. 
The solution is unique in either case. 
Proof. 

(a.) The solution exists, with the indicated properties, because the equation admits a subsolution (^_, with < if-, 
which is the solution of the equation (from section 7) 

A^ip^ - r(7)(p_ - bip'!^ = 0. (83) 

The solution satisfies (f- < 1 under the hypothesis made on r{'~f) because the equation for ip^ admits then a su- 
persolution equal to 1. The original Lichnerowicz equation (81) admits as supersolution ip^ > 1 the solution of the 
equation (cf. section 8) 

A^ip+ + aip:^:^ + q^^^' = (84) 

because we have 

r{j)p+ +bip'!f_>0 if ip+>l and r{j) +b>0. (85) 

(b.) When (M, 7) is in the positive Yamabe class, the equation is equivalent to an equation of the same type with 
zero linear term because of conformal covariance. We may then argue existence just as in (a.), because the condition 
on r(7) when it is negative has become vacuous. 

The solution tending to 1 at infinity of the equation with r{'j) = is unique because of the monotonicity of / in 



if. In the general case one uses the conformal transformation of curvature as in Section VI . Take for simplicity of 
writing q=0. We have now if tpi, z = 1 or 2, is a solution, 

-r{^?l) = -b + a^rP-Q; (86) 

therefore the conformal identity with u = ip'^ p2 gives 

A^2p^M +{b- a</77^"^+^V =ib- aip~^^'^^^)u^ . (87) 

This equation may be written 

^^r." - 1" (^^) + ■'^'"'*'' (^:^) }"<■"- »-o m 

If u > 0, 6 > 0, a > the function u, which tends to 1 at infinity, can only be u = 1 on M. 

Remark 1. We see that the condition that (M, 7) be in the positive Yamabe class is not necessary for the existence 
of a positive solution if 6 ^ 0. However if 6 ^ the Hamiltonian constraint is coupled with the momentum constraint, 
and its solution is not the whole story. 

Remark 2. The condition b > — ^(7) will somewhat be relaxed in the last section but we will require b > 0. 

X. DISCONTINUOUS SOURCES 

It is essential for physical applications to admit isolated sources, hence discontinuous functions q. This possibility 
is included if we extend the previous existence theorem to functions q e VFq 5+2- We also will take a € W^g^2 to 
include the possibility of discontinuous scaled momentum v. We take d = b ^ ^s s+2' ■s > "■ We leave more general 
cases for later study. 

Theorem. The Lichnerowicz equation with scaled sources, 

-p' 



A^ip ~ r(7)<^ = /(., ip) = -a<^-^ - q^-^ + bp'^ , (89) 

on a, {p, a, p), a > - + 2 , p > —- asymptotically euclidean manifold (M, 7) in the positive Yamabe class has one 
and only one solution ip > , p — 1 = u e W^ ^ if a, 9 G ^05+2 '^i*^ 5 > —- , p > ^ , a>0, g>0, 6>0, 5e 

Proof. We first conformally transform the equation to an equation with no linear term. We then proceed as follows. 
Consider a Cauchy sequence a^, q^ £ ^^5+2' s > —, converging in the VFj^- , 2 norm to a, q. Denote \>y (p^ = \ + u^ 
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the solution with coefficients a^, q„. We know that u^ G Wj_^_2 g and that there exists numbers £ > (depending only 
on (M, 7) and b) and m> £ (depending only on {M, 7) and the Wj^ , 2 norms of a and q) such that £ < Lp^ <m. 
The difference Un — u^ satisfies the equation: 

A^(u^ - u^) - yl^^(u^ - u^) = ^'^^{ai, - a^) + ipj^^ {q„ - q^) (90) 

with: 

/ ,n-P - in-P\ / ,n-P' - ,n-P'\ /,nQ-,nQ\ 

(91) 

Recall that for n = 3 we have P = 7 , P' = 3 and Q = 5. The quotients in the above formulas are then polynomials 
(with coefficients equal to 1) in (pT^^ and (/?^^ for the first two, and (^^ and (p^^ for the third. Therefore, they are on the 
one hand positive and, on the other hand, uniformly bounded (for any pair v, 11) because < £ < (pfj_, (fi, < m. For 
general n the numbers P, P' and Q arc positive rationals, the quotients in the formula are also positive and uniformly 
bounded. We deduce from this uniform boundcdncss that there exists a number N such that 




\A 



^-llw-oV. ^ ^{ll«-lll^oV2 + ll'^-ll<. + 2 + ll^ll<,. + j • (92) 



We infer from this estimate and the positivity of A^i, that the operator A^ — A^^ is injective in W2 g (See Condition 
2 in Theorem 1 of Appendix B). Therefore, there exists a number C depending only on (M, 7), the VFj^ 1 2 norm of 
6, and the VFj^- , 2 norms of a and q such that 

Iki. - Uf,\\wl, < CJllaj, - af^WwS,^^ + \K - (If^Ww^l,^^] ■ (93) 

Since ai, and q^, are Cauchy sequences, the same is true of u^, because of the above inequality. Hence Ui, converges 
in W2 g to a limit u G W^f^- The convergence is a fortiori in C" if p > ^, and for some positive a, since S > — — . 
The u^'s are such that I + u^ > £ > 0; therefore also (p = l+u>£>0. The function ip satisfies the Lichnerowicz 
equation (in the sense of generalized derivatives) with scaled sources. 

XI. GENERAL CASES 

In the case where there are unsealed sources the coefficient d in the Lichnerowicz equation is negative or zero on 
a maximal initial manifold M. It can take different signs if M is not maximal. The previous simple method to 
obtain sub and super solutions does not apply. We will then look for constant sub and super solutions £ and to, 
< ^ < 1 < TO. We will also obtain a new theorem for the Lichnerowicz equation in the case of scaled sources on a 
non maximal submanifold. To make the algebra easier we restrict our study to the important physical case n — 3. 
Results along the same lines can likely be obtained for general n. The equation is then 

A^p — rip + ap> + qp}^ — dp — 0, 

a > 0, g > 0, d^b-c, 6 > 0, c > 0. (94) 

The numbers £ and to are admissible sub and supersolutions if they satisfy on M the following inequalities: 

Pxi£'^)<0, Pxim'^)>0, for aUxe M, <£< 1 < TO (95) 

where Px is the polynomial 

Px{z) = d{x)z^ + r{x)z^ - q{x)z - a{x). (96) 

Remark. In the case of n > 3 the problem is the study of the sign of the function: 

Fx{z) = d{x)z'^ + r(a;)z"-i - q{x)z^''-^^''^ - a{x) (97) 

for numbers £4/("~2) ^^j^^j j^4/(n-2) 

Since all the coefficients in Px tend to zero at infinity the conditions that we will obtain depend on the ratios of 
their respective decays. 
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We denote by M+ the subset of M where d > 0, by M_ the subset where d < 0, by Mq the subset where d = 0. 
In the case of isolated sources -M_ is a compact subset of M . We study the sign of P^ on these various subsets. The 
derivative of Px is 

dP^/dz ^ M{x)z^ + 2r{x)z - q{x). (98) 

1. On M+, d{x) > 0, the derivative dPx/dz has 2 roots of opposite signs. The positive root is 

c,(.) ^ ^^:M±m±^^^ (99) 

We have Q+{x) > if r{x) < 0, or if r(x) > and q{x) > 0. 

dPx/dz is equal to —q{x) < for z = and is negative or zero as long as z < (.^(x). Therefore P^ decreases from 
a{x) < for z = to a minimum for z — C+l^^) ^md then increases to +oo when z increases to +oo. Hence Px has 
one and only one positive root z+{x). We have Px{z) > as long as z > z+(x). 

There exists £{x) > such that Px{£^{x)) < if and only if z+(x) > 0. Indeed numbers i{x) and m{x) such that 

< £{x) < z+{x) < m{x) , xe M+ (100) 

satisfy 

Pxii^ix)) < 0, Px{m\x)) > 0. (101) 

Lemma 1. There exist numbers £+ and to+ such that: 

PxiiX) ^ 0' Pxim\) > for all x e M+ (102) 

if and only if 



and 



inf z+(x) > 0. (103) 

xeM+ 



sup z+(x) < +00. (104) 

xGM-i- 



^iimum^m"']>o im 



Sufficient conditions for the first inequality are 

inf JZlM 
x^iii+] 3d{x) ' y9d'^{x) ' 3d{x) 



inf ---^M-->0. (106) 

xeM+ d{x) + \r{x)\ 

For the second inequality they are that -^S^, fS'j ^t^t be uniformly bounded on M+. 

Proof. The necessary condition as well as the first sufficient condition are consequences of the previous study. 
Sufficient conditions for this first condition to be satisfied are that one of the two terms in the sum has a strictly 
positive infimum. The second sufficient condition results from the fact (elementary calculus) that Px{z) < if 

z<minfl "["•^ ) (107) 

V d{x) + \r{x)\J 

Remark. The sufficient conditions will be satisfied on the whole of M+ if we can split it into two subsets, M_|_ = 
Ml U M2, such that 

-f 44^^ > and inf ^ > 0. (108) 

xeMi d{x) + \r{x)\ xeM^ d{x) 
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This pair of inequalities can be realized when M is compact and a{x) + q{x) ^ by a conformal change of choice of 
the metric 7 to a metric 7' having a strictly negative curvature in the complement of Mi in M . Such a construction 
can also eventually be made in the asymptotically flat case, by resolution of an adequate Dirichlet problem. 

2. On M_, d{x) < 0. 

We have Px{z) < for all z > 0, hence no admissible m{x), if r{x) < 0. We therefore suppose r{x) > for all 
X e M_. If r^(x) + 'Sq{x)d{x) < 0, we have dPx/dz < for all z; and the polynomial P^ takes non negative values 
only if it is identically zero. If r^(x) + 3q{x)d{x) > the polynomial dP/dz has two positive roots: 

, ^ r{x)-{r^{x)-3q{x)\d{x)\y^^ , ^ 

3|rf(a;)| 

C2X = ^ ^ ^ "- I y ^' ^ ^'^ >0, 110 

with Ci(a;) > if and only if q{x) 7^ 0. 

The polynomial P^ decreases for < z < (,\{x)^ increases for (,\(x) < z < C,2{x), and decreases to —00 for z > C2(2;). 
We have Px{0) ~ —a{x) < 0. Therefore P^ takes negative values for some z > if either a{x) > or Ci(a;) > 0, i.e. 
q{x) > 0. The polynomial P^ takes positive values, equivalently admits two positive roots zi(x) and 22(2;) which are 
such that 

Cl{x)<Zi{x)<C2{x)<Z2{x), (HI) 

if and only if its maximum, attained for z = C2(a;), is positive, 

Px{C2{x)) > 0. (112) 

We have then Px{z) < for < z < zi(x), and Px{z) > for z > Z2{x). If 

r^{x) + 3q{x)d{x)<0, (113) 

the polynomial Px is always decreasing. It takes positive (i.e., non-negative) values only if it is identically zero. 

Lemma 2. Suppose that r{x) > 0, r^(x) — 3q{x)d{x) > and Px{C2{x)) > for all x £ M^. There exist numbers 
i^ and TO_ such that: 

Px{ei)<0, Px{mi)>0 forallxeM_ (114) 

if the following conditions are satisfied: 

inf zi(x) > 0, sup zi(x) < inf Z2(x), (115) 

M_ xeM- x<EM_ 

and ^, ^, ^ are uniformly bounded on M_. 
Proof. All numbers i- and ?Ti_ such that: 

e^<zi{x), zi{x) < m- < Z2{x) ioT all X e M- (116) 

are such that Px{f-^) < 0, Px{m^) > 0. These numbers exist, with £_ > and +00 > m_ > £- under the given 
conditions. 

3. On Mq, d{x) — 0, Px reduces to a second order polynomial 

Px{z) ~ r{x)z — q{x)z — a{x). (117) 

If r{x) < 0, then P^ < as soon as z > except if it is identically zero. We suppose r{x) > 0. Then Px admits one 
positive root zo{x): 

zoix) = (2r(x))-i {q^{x) + 4a(a;)r(a;)}'^^ > 0. (118) 

Lemma 3. We suppose that r{x) > for all x e Mq. 

There exist £a > and mo > £a such that Pxi^o) — and Px{TnQ) > for all x E Mo if and only if 
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and 



inf 44 > or inf 44 > (119) 

x^Mq r{x) x£Mq r{x) 



Gl X ) oix] 

sup —-^ < +00 and sup — ^— < +oo. (120) 



Proof. Under one or the other of the first inequahties we have 



The other ones insure 



We set 



inf zo(x) > 0. (121) 

xeMo 



sup 2:0(2^) < +00. (122) 

xGMo 



^Q = inf zq (2;) , TOq = supzo(a^)- (123) 

Mo Mo 

All numbers i and m satisfying the following inequalities 

< ^ < ^0 = inf zo(a;) < sup 2:0(2;) = mo < m (124) 

Mo xeMo 

satisfy Px{£'^) < and Px{m'^) > for all x G Mq. 

The following lemma is an immediate consequence of the previous three. 

Lemma 4. We suppose that the conditions given in the lemmas 1, 2, and 3 for the existence of £_, i+, £0 and 
TO_, 771+, 777o are satisfied. Then there exists £ and m such that: 

<^< 1 <7n and P:,(^'*) <0 , Pxim^) > for all x e M (125) 

if the following inequalities hold: 

771+ < 77l_ , Too < 77l_ , 771_ > 1 . (126) 

Proof. Take 

m = m^, £ = min(l, 4, £+,£_). (127) 

Then £ and to, satisfy the required inequalities for all x g M . They are admissible sub and supersolutions. 
Theorem. On a 3-dimensional asymptotically euclidean manifold the Lichnerowicz equation 

A.yip — rip + aip^'' +q(f — d(f = , 
a>0, g>0, d ^ h-c, 5>0, c>0. (128) 

with r, a, q, d £ Wfg^^^ s > -, ~- < 6 < n — 2 — -, admits a solution (p > 0, (/? — 1 G Wf+2 s ^^ *^^ assumptions of 
Lemma 4 are satisfied. 

Corollary No Unsealed Sourees, d=b>0. The Lichnerowicz equation has a solution ip > 0, ip — 1 ^ ^s+2 s i^ 

(i) The quotients -jr^, frfy, fT^y are uniformly bounded. 

(ii) There is a positive number e > such that if l(J\ < e, then 

rix) < and ^T^^ > e' > 0. (129) 

b{x) 

This last condition can be achieved if a + g ^ by a conformal transformation and solution of a Dirichlet problem in 
the subset of M where (a + q)/d < e, so long as this subset is compact (cf. [|]J|]). 
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XII. UNSCALED SOURCES, CASE A^ = 3 

We treat in this section the haniihonian constraint for unsealed sources in the case n — 3. The Lichnerowicz 
equation reads 

A^f-rip + aip^'^ + cf'^^0 (130) 

The functions a > and c > are given on (M, 7). 

Theorem. Let (M, 7) be a (p,(T,p) asymptoticaUy eucUdean manifold, a > - + 2, a > —- with r > 0. Let 
a, c £ W^ g,2 be given on (M, 7), s > —, a > — — . There exists an open set of values of a and c such that the 

Lichnerowicz equation with unsealed sources has a solution (f > 0, with 1 — ip E T^f+2 s- 
Proof, We look for constant admissible sub and supersolutions i and m such that 

<£<l<m , 
Pxii^) > 0, Pxim^) < 0, for all x e M, (131) 

where P is the polynomial, 

P^{z) = c{x)z^ -r{x)z^ +a{x). (132) 

1. Case c > 0. 

We set z ^ X^^ and consider the polynomial which has the same sign as P^, 

Q{X) = a{X^ -a-\X + a-^c}. (133) 

This polynomial has 3 real roots if 

4r^ > 27ac^ . (134) 

Two of these roots are positive, given by the classical formulas: 

e ^ e + 2n , , 

X2 — A sin- , Xi — A sin — - — (135) 



with: 



-(S^' --^^^ »^«4- 



The corresponding roots of P^ are 



zi = ^<^ = Z2ix). (137) 



Xl - X; 
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We have Px{z) > for < z < zi(x), P^{z) < for zi{x) < z < Z2{x). 
2. Case c{x) = 0. 
The polynomial Px reduces to: 

Px{z) = -r{x)z^ + a{x). (138) 

We have Px{z) > for < z < {r^'^aY/'^, and P^iz) < for z > (r^'^aY/'^. Note that (r^ia)^/^ jg ^^ve value for c = 
of the previously computed zi while the previous Z2 tends to infinity when c tends to zero. The eases c{x) > are 
thus unified. 

The following constants £ and m are sub and supersolutions if 

£<zi{x), zi{x) < m < Z2ix) foraUxeM. (139) 

They exist, satisfying the required properties < ^ < 1 < ?7i, if 

inf Zi(x) > , inf Z2(x) > max < 1, sup zi(x) > . (140) 

x&M xeM \ ^gM J 
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One then takes 

i — niin < 1, inf zi (x) > , m — inf Z2(x). (141) 

[ xeM J xeM 

We give below sufficient conditions to satisfy the various inequalities, using the expressions 

2 V a 2 V CL 7T 

Zl = ^ . ,,^ ^r^TTTT, Z2 = ^^^TTTTrr, withO<^<-. (142) 

V3sin((6' + 27r)/3)' V3sin(6l/3)' " "2 "■ ' 

The functions sin(0/3) and sin((6' + 27r)/3) are respectively increasing and decreasing when 9 increases from to ^. 
Denote by ^min and ^max the infimum and supremum of on M defined as solutions between and ^ of the equations 

. ^ . ^3c(3a)i/2 . 3c(3a)i/2 

sm 6^ini„ = mf — - — -— — , sm 6^,nax = sup — -p—. (143) 

M 2rr^i^ M 2rr^/^ 



Therefore: 



inf 01 > ^ inf (r-ia) {sin((6'^in + 27r)/3)}"^ 
^ v3 ^^ 

2 
supzi < — =sup(r"^a) {sin((6'max + 27r)/3)}"^ 

M V3 M 

inf Z2 > 4Tii^f('^"'«) {sin^max}"' . (144) 



We find by elementary calculus that: 



< sin((0„,ax + 27r)/3) - sin(0„iax/3) = \/3cos((0„,ax + 7r)/3) < ^ (145) 



3c(3a)i/2 
^^P o 1/2 = 1 • (146) 

To insure the existence of constants i and to satisfying the required properties we suppose 

1^,1.... 

V3 M 



mf(r-ia){sin((0,„i„ + 27r)/3)} '>0, (147) 



and we set 



We suppose also 



and 



£ = min{^_,l}. (148) 



m+ = ^inf(r-ia){sin(0„ax/3)} '>1 (149) 



V3 M 



m- = —= sup(r-ia) {sin((6l„,ax + 27r)/3)} ^ < to+ . (150) 

V3 M 

The condition ?ti_ < to+ can be satisfied for an open set of values of the coefficients c, a (given r) due to the previous 
elementary study. We can take for m any number between max{l,TO_} and to,4_. The numbers £ and m so chosen 
are admissible sub and supersolutions of the Lichnerowicz equation. The existence of a solution ip with the required 
properties results from the general theorem, given in Section 5. 
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XIII. COUPLED SYSTEM 

In the conformal method the momentum and the hamihonian constraints decouple when the initial manifold M has 
constant mean extrinsic curvature and the unsealed sources have a momentum A'^ = 0. The theorems of the previous 
sections are then sufficient to give existence, non-existence or uniqueness theorems of the systems of constraints. 
We will in the next sections study cases where one of these hypothesis does not hold; hence the constraints do not 
decouple. 

XIV. IMPLICIT FUNCTION THEOREM METHOD 

The use of the implicit function theorem is the simplest way of proving existence of solutions of equations in the 
neighbourhood of a given one. It works as follows. 

Let U and V be open sets of Banach spaces X and Y and let J- he a, C^ mapping from U xV into another Banach 
space Z: 

J^:UxV^Zhyix,y)^J^{x,y). (151) 

Suppose that the partial derivative of J^ with respect to y at a point (xq, yo) lE U x V, T'Axq, j/q), is an isomorphism 
from Y onto Z\ then there exists a neighbourhood W of xq in U such that the equation 

.F(x,2/)=0 (152) 

has a solution y e 1/ for each x E W. 

We consider the quantities q and v (scaled sources) together with N and u, a traceless symmetric 2-tensor as given 
on the asymptotically euclidean manifold (M, 7), with q, v, 1 — N~^ G ^sS+2' "^ ^ ^f+i i5+i- ^^ ^^^^ discuss the 
existence of ip and /3 as we perturb J and t away from zero. The points x, y and the Banach spaces X, Y, and Z are 
as follows: 

x^{T,j)eX^w:^,^,^,xWl,^,, 

2/ ^ (/3,^- 1) e ^ ^ y n {^ > 0}, y ^ w^^.^s x w!+2,s ■ 

Z ^ Wl,^^ X Wls^, . (153) 

The mapping J-' is given by 

J^ix, y) = (H(t, J; ^, /3), M(t, J; ^, /3)) (154) 

where TL and AA are the left hand sides of the conformal formulation of the constraints, T-i(x,y) = A^ip — 
/(.; r; /3, p), M{x, y) = V.{{2N)-^CI3) - h{- r, J; p). 

The multiplication properties of weighted Sobolev spaces show that JF is a C^ mapping from X xV into Z ii s > - 
and S > — — . The partial derivative J^' at a point (x,y) is the linear mapping from Y into Z given by: 

{Sp,Sp)^{n'y,M'y).{S(3,Sp) (155) 

with {A is given by ([iq)) 

n'y.iSP, 6ip) = A^dip - aSip + 2knP^^2N-^A.CSp , 

a = r + Pp-P-^a{(3) + P'qp-^'-^ + dQp^-^ , 

M'y.{5p,5p) = \/.{2N-^C5(i) - \5p , 

(n — 2) n — 2 



(156) 



Theorem. Specify on the asymptotically euclidean manifold (M, 7) a traceless tensor u £ W^f 1 1 511, a scalar A^ > 
with A^ — 1 e W^f+2 5; and scaled and unsealed sources q, v, c G W^g,2i s > -, —— < S < n — 2— -. Let {/3q, pq) be 
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a solution of the corresponding constraints with Dtq = (hence tq = since tq G T^f+i 5+1)1 ^^^d Jq = 0. Suppose 
that on M 

ao = r + P^^P-^a{po) + PV(7^'"'(Z - cQ^'i'^ > . (157) 

Then there exists a neighbourhood U of (tq, Jq) in X such that the coupled constraints have one and only one solution 
{(3,^), ip>0, {(3,l-^)eY. 

Proof Under the hypotheses that we have made the partial derivative J^y{xo,yo) is an isomorphism from Y onto Z 
because the system of linear elliptic equations 

V.{(2iV)-i£(5/3} = h 

A^S(p - aaSif = -2knip-^{2N-^A).C5(3 + k (158) 

has one and only one solution {5ip, 5j3) G Y for any pair {h, k) G Z . 

Corollary. The conclusion of the theorem holds if (M, 7) is in the positive Yamabe class and d > (realized in 
particular if all sources are scaled) without having to consider the sign of a^. 

Proof. If (M, 7) is in the positive Yamabe class we can choose (M, 7') in the same conformal class and such that 
r(7') — 0. To the data N, u, q, v correspond data N' , u', q' , v' and to the solution /3o, (fo corresponds a solution of 
the transformed conformal constraints. The corresponding Uq is positive and the conclusion of the theorem applies 
to the transformed system, hence also to the original system. 

XV. CONSTRUCTIVE METHOD WITH SCALED SOURCES 

We will give in the next two sections another method to obtain solutions of the coupled system. It will give new 
results for unsealed sources on a maximal manifold. It is possible, though not proven, that the hypotheses we make 
in the case of scaled sources on a non-maximal manifold imply that this manifold is in the positive Yamabe class. 

Lemma 1. The equation 

A^Lp — f{.,ip) = rep — aif^ — qip^ + bip^ (159) 

with r, a, q, h satisfying the hypothesis of the theorem in Section 9 admits as a supersolution the solution $(A), 1 — 
^{A) G VKf+2 S-i of t^'2 equation 

A^^^ fA{;^)^-A^-P -q^-P' (160) 

if a < A, with A a given function in Wj^ , 3. 

Proof. The function ^{A) exists by the theorem in Section 9. It satisfies 

A^$ - /(., $) - /a(., $) - /(., $) < (a - A)$-^ < , (161) 

hence it is a supersolution. 

Theorem. Under the conditions on r and h given in the theorem of Section 9 there exists a number e > such 
that if 

\\Dt\\wI,^^ < e , n>p,5> -n/p, (162) 

the coupled constraints admit a solution (/3, ip) with /3, \ — ip £ W^f+2 s- 
Proof. We will construct a sequence {(p^,f3„) by the inductive algorithm 

A^y^i. = f {.,13^-1, p^) = ript, - a(/3i._i)(^^^ - q'-PZ^ + ^V'? , 



V^{{2N-^){Lp,r] - h\.,^,) = V, {{2N)-^u'=] 



n 
with 
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+ 1^2„/(«-2)yj^ ^ ^J (^g3) 



n — 2 

a{l3) = k,,{2N)-^\-u + CpW fc„ = - -. (164) 

4(n — i) 

The equations for the (p^s admit all the same subsolution (/?_, which depends only on r and h. They admit the 
same supersolution $(^) if there exists A e W^ ^^^ ^^^h that a{(3v-i) < A for all u. 
We start for instance from /3o = and choose A such that 

A > a{Q) = kn{2N)-'^\u\^ . (165) 

Suppose that fiu~i G Wgg^2 ^^'^ that a(/3^_i) < A. Then (pi, exists, 1 — ip„ e ^s+2 5' ^^'^ P- < Pv < ^(^)- Also, 
Pu £ W^,2 s exists, and there is a constant C (cf. Appendix A) depending only on (M, 7) such that 

II/3.||h.|,, <C^{;^ll^r||^„^,,,sup$(A)-/("-) + ||P||^.^^^} , (166) 

where P is the given vector 

P = ^.{{2N)-^u}+v . (167) 

The weighted Sobolev multiphcation theorem and the expression for a{(3) imply that if /3^ e VF^^^, ^ < 1, 5 > — ^, 
then 

a(f3^) e VFf „ , „ for aU 6' such that 5' < (5 + (5 + -) , (168) 

and there exists a number C such that 

||a(/3.)|lw. ,^^ < C|||^i||i. ^ +I1/3.I1Lp 1. (169) 

1,0+2 ^_ 1,0 + 1 2, a J 

By the weighted Sobolev inclusion theorem, there exists then another constant C such that 

||a(/3.)|lco„<C||a(/3.)||i^.^,^^ (170) 

for aU S" <S' + 2+^, hence also for aU S" <2S + 2 + ^ . 
Since S + — > there exists a number a such that 

Tl Tl 

6 + 2 + -<a<2S + 2 + -. (171) 

P P 

We choose for A a function of the form, with fi some positive constant, 

A = ^l/a" , (172) 

where a = 1 + (P (see section 3). We have A £ W^o 5+2 since a > S + 2 + - . For such a function A, the inequality 
fl(/3i/) < v4 is equivalent to 

a"a(/3,)<^, i.e., |la(/3,)||^2^„<^. (173) 

Using the previous estimates we see that a sufficient condition to insure on M the inequality a{l3^) < ^ is to have 
some number depending only on (M, 7) and A^, denoted by C, that has the property 

where we have set $^ = $((T^"/i). 
We choose /x large enough to have: 

C^i>S, S=\\v\\l,. +\\u\\l,. . (175) 

The inequality obtained above shows that we can construct (pi,+i and hence (3^+i, enjoying the same properties as 
(fiy, Pi, if Dt is sufficiently small in VF^^- 1 2 norm. The existence of a solution (<p,/3) of the coupled constraints as 
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limit of a subsequence is proved by a compactness argument and elliptic regularity as in the case of the hamiltonian 
constraint with a given [3. 

Remark. The number e depends on the choice of /i and the function $^. Neither of those depends on r or on 
h, i.e., on t. However, our theorem imposes a restriction on the size of t on the subset of the manifold M where 
r < 0, since it supposes that b = -^^^t^ > —r. This could lead to a difficulty, pointed out by N. O'Murchadha, on an 
asymptotically euclidean manifold where r would be too negative. Indeed it is known that in the case p = 2, S = —1 
there exists a constant Cp such that the following Poincare estimate gives an upper bound of |t| in terms of \Dt\: 

\\r\\Ho.^,<Cp\\Dr\\H,,„. (176) 

If we suppose that an analogous inequality holds for p > n, S > — — , i.e, that there exists a constant Cp such that 

\\r\\^^^^ < Cp\\Dr\\^.^^^^ , (177) 

then by using the Sobolev embedding theorem: 

supCT^^r^ < CCI\\Dt\\1,p , for all a < 5 + - . (178) 

j^f f" "''^0.6+1 p 

This inequality implies that the condition b > — r can be satisfied only if r satisfies the condition 

- ra2" < -^CC^e . (179) 

71 — 2 

We can estimate the value of e as follows, considering for simplicity the vacuum case q — — v — 0. The supersolution 
<I>^ satisfies the equation 

A^<i>^ = ~a<d;^ (180) 

that $,, > 1; therefore A^^^ < A ; 
equation 



We know that <i>^, > 1; therefore A^^^ < A and $^ < vp^, where ^^ is the solution with *,, - 1 G W^^2 s °f ^he 



A^^^ = -A = -m/ct" . (181) 

Obviously ^^ = 1 + ^iWl, where wi depends only on (M, 7), satisfies the equation 

A^wi = -1/cr" (182) 

and tends to zero at infinity. The inequality to satisfy is then 

WDtWLp < (C/i-S')(l + Ci/i)"'*"/^""^\ withCi^supwi . (183) 

0,S + 2 ^j 

The right hand side is maximum for a finite value ^ — fiQ with fiQ given by 

in~2)C + AnC^S 
^" = {3n + 2)CC, ■ (^''^ 

(Note that fiQ > -^ ii n > 2.) We find therefore 

4„ A.a5M-'"-^'"-2^, cy j^^^, 



C J j 3n + 2\ Ci 

It is an open problem to prove the analogue of the Poincare inequality in W^ g and to decide whether the restriction 
imposed on r implies that (Af, 7) is in the positive Yamabe class or not. The conclusion may be (cf. H]) related to 
the existence of apparent horizons as in the proof by Schoen and Yau M] of the positive energy theorem. 
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XVI. COUPLED CONSTRAINTS WITH UNSCALED 
SOURCES 

We treat in this section the system of constraints for unsealed sources on a maximal submanifold in the case n = 3. 
This system is coupled if the given initial momentum J does not vanish. The equations are: 

A^ip — rip + a{l3)<p)^'^ + ap^ — 
\I^{{2N)-^{CPY^] = Vj{(2iV)-iM'J}+^iojJ'. (186) 

The functions c > and A^ > and the tensor u and vector J are given on (M, 7) . We denote by /3o the solution of 
the equation: 

y,{{2N)-\CI3^r} = ^,{{2N)-\'^} . (187) 



We have the following straightforward result. 

Theorem. We suppose that the given quantities r, c, J € W^ ^^2^ " '= ^f+i s+n I — N e W^^2 s^ P > n, S > 

s,S+2 



—n/p are such that there exist positive functions A^, A^ G ^sS+2 with: 



A_ < a(/3o) < A+ (188) 

and such that the equation: 

A^(p -rip + A+ip-^ + ap^ ^ (189) 

admits as supersolution the constant rn^^ > 1, and that the analogous equation constructed with A_ admits as 
subsolution the constant £a_ < 1, ^a_ > 0. Then there exists e > such that 

\\J\\ws,,,<^ (190) 

implies that the coupled constraint equations have a solution (/3, tp) with (/? > and /3, 1 — (p € W^,2 5- 

Proof. It is elementary to check that iA_ and niA^ are admissible sub and supersolutions of the Lichnerowicz 
equation constructed with a{(3) if 

A- < a{(3) < A+ . (191) 

In this case the equation has a solution ip with iA_ < ^ < rriA^ ■ 

The momentum current being independent of other quantities its bound does not affect other estimates. 

We will construct a sequence ^PutPu as in the previous section. We now have to show A- < a{(3,^-i) < A+ implies 
the same inequalities for a(/3i,) if (/3i/_i < rnA+ and || J||vf'' is small enough. We use the fact that (the notation | | 
means here the 7 norm): 

|{a(/3o)}i/2-(167V)-i|/:(/3-/3o)|| 

< {a{l3)Y'^ < {a(/3o)}i/2 + (16iV)-i|£(/3 - /3o)| . (192) 

We deduce from the momentum constraint satisfied by /3y the elliptic estimate 

11/3. - Mwl, < Cm]^JJ\\wl,^, • (193) 

The proof is then completed along the same lines of previous proofs. 
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XVIII. APPENDIX A: ELLIPTIC LINEAR SYSTEMS ON MANIFOLDS EUCLIDEAN AT INFINITY 

For the convenience of the reader we recapitulate here some known facts. 

A hnear differential operator of order m from sections u of a tensor bundle E over a smooth Riemannian manifold 
(M, 7) into sections of another such bundle F reads 



Lu = J2'^kD''u (194) 



fc=0 

with o/c a linear map from tensor fields to tensor fields given also by tensor fields over M. 

The principal symbol of the operator _L at a point x S M, for a covector ^ at x, is the linear map from E^ to F^ 
determined by the contraction of a^ with ((8)0™. The operator is said to be elliptic if for each x G M and ^ G T*M 
its principal symbol is an isomorphism from E^ onto F^ for all ^ 7^ 0. 

Example. The conformal Laplace operator in a metric 7 on M acting from vector fields /3 into vector fields is 



W,{Cpy^ = V, (w'f3^ + V^P' - '^f^Vkl3'''\ 



(195) 



Its principal symbol at x, with ^ e T^M , is the linear mapping from covariant vectors b into covariant vectors a given 

by 

e^^b,+e^A-l^l^'bk^a,. (196) 

This linear mapping is an isomorphism if ^ 7^ because 

2 



{Cmb'b,)+[l--j{Chr>0 (197) 

if ^ 7^ and 6 7^ 0. The conformal Killing operator is elliptic. 

Theorem. Let (M, e) be a (complete) Riemannian manifold Euclidean at infinity. Let 

m 

Lu = ^akD^u (198) 

fc=0 

be an elliptic operator on (M, e) . Suppose the coefficients of L satisfy the following hypotheses 

1. There is a C°° tensor field Am on M, constant in each end of (M, e) such that for some p with 1 < p < +00 



n _ ^ n 



am- Am ^yy^ g , Sm > - + 1 , Om > , 



a-k S W^ : ,Sfc> — hfc = m + l,(5fc>TO — fc ,0<fc. 



2. 



Then for each s such that Sk + m > s > m the operator L maps W^ § into VFj_„j ^ 1 ,„ with finite dimensional kernel 
and closed range if 

Ti 71 

<d < -m + n . (199) 

P P 

If, moreover, L is injectivc on W^ g then it is an isomorphism and there is a number C such that for each u in W^g 
the following inequality holds: 

llwll^P^ < C\\Lu\\^._^^^^^ . (200) 

This theorem applies to the Poisson operator A — A: under the hypothesis indicated in the theorem in Appendix B. 
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XIX. APPENDIX B: SOLUTION OF A^ip = F{X,ip) ON AN ASYMPTOTICALLY EUCLIDEAN {M,j) 

Let A-y denote the Laplace operator on scalar functions on (M, 7). Let / be a real valued function on M x /, with 
/ an interval of R, given by {x,y) h^ f{x,y). We will show that the sub and supersolution method used by one of 
us (J.L) in the case of a compact manifold can be extended to asymptotically euclidean ones. Recall that (M, 7) is 
a (p, cr, p) asymptotically euclidean manifold M of dimension n, if 7 — e G W^ ^ with (M, e) euclidean at infinity and 

1. Linear Equations 

Definition. Suppose M is not compact. We say that / tends to a value c G M at infinity if for any e > there exists 
a compact 5* such that: 

sup 1/ - c| < e . (201) 

M-S 

Lemma 1 (maximum principle). Let (M, 7) be an asymptotically euclidean manifold. Suppose that a C^ function 
if on M satisfies an inequality of the form: 

A^ip + a.Dip-hip <0 (202) 

with a dot denoting the scalar product in the metric 7, while a and h are respectively a vector field and a function 
on M, both bounded. Suppose that h> on M . Then 

a. If (p tends to c > at infinity then there exists a number £ > such that (p> i on M . 

b. If (p tends to c = at infinity then </? > on M . 

Proof. One knows by the classical maximum principle that if ip attains a nonpositive minimum A at a point of M 
then (/3 = A on M . Also, if Z) is a bounded domain of M with smooth boundary dD and if the function p attains a 
nonpositive minimum in D U dD this minimum must be attained on the boundary dD. 

a. Choose e > so small that e < c. li p tends to c at infinity there is a compact S such that p>c~e>Qon 
M — K. Imbed 5 in a relatively compact domain D with smooth boundary dD. On dD, p takes positive values, 
therefore p does not attain a nonpositive minimum on the compact set D U dD; it attains a positive minimum c'. 
The number d is the smaller of the two positive numbers c — e and c'. 

b. Suppose that p takes a negative value a on M. Choose e < |a|. There is a compact S such that 

sup I/I < e . (203) 

M-S 

Take a relatively compact open set D containing 5*. If p takes a nonpositive minimum it is on the boundary dD, i.e., 
in M — S, which contradicts the fact that the absolute value of this minimum is necessarily greater than or equal to 
|a|, itself greater than e, which is the maximum of \p\ in M — S. 
Theorem. Let (M, 7) be an {p,a,p) asymptotically euclidean manifold. Let k E W^g,2J S > —- he given. The 



operator A-y — fc is injective on Wji 2 § if either 



1. fc > , s > - . (204) 

P 

2. / {\Df\^ + kf}p, > for all / e C^ , / ^ , 
Jm 



IM 

Ti ft 

s>0and(5> I Up ^2, 5>-liip = 2 

2 p 



(205) 



Corollary. Under the hypotheses 1 or 2 the operator A-y — fc is an isomorphism from Wj, 2 5 onto W^ g,^ if 

s < cr - 1 - 2:, -2 < (5 < -'s- + n - 2. 

— p ^ p p 

Proof. 1. If s > -, a solution in VFf^_2 g with S > —- is in C^ for some positive a. The difference 7 — e is in Ci for 
some positive p. The maximum principle applies and shows that u = on Af. 

2. The solution u € W2 g is not necessarily C^. To prove that u = wc will multiply by u the equation and 
integrate on M. 

If mGC^, then 

uA^up^ — — Du.Dup,^ . (206) 

M Jm 
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We can estimate the integrals in the above formula in terms of the Sobolev norm. In the case p = 2 we have 

uA^ufi^ < ||7/||ff„ JjA^u||H„,^,sup(l + d2^-(*+i) , (207) 

M M 

which is a bounded quantity whenever 6 + 1 > 0. 
In the case p y^ 2 we have 



The considered L^ norm is bounded if 

'r^ -t- 1 "i > r?, i p A + 1 ^ 

2 p 



Tl Ti 

2p'{S + l)>n, i.e. , 5+l> -~- . (209) 



The same kind of estimate applies to the integral of Du.Du. 

The density of C^ in the weighted Sobolev spaces shows then that a solution u G W^ g satisfies the equality 



uA.yUfi~f — — Du.Dufi~f — / ku fj,.y . (210) 

M JM J M 

Therefore, under the hypothesis made, we have Du = 0; hence u =constant and u — Q because u tends to zero at 
infinity. 

The corollary is a consequence of the general theorem on elliptic systems on an asymptotically euclidean manifold 
recalled in Appendix A. 

If n > 2 the inequality (5<— - + n — 2is compatible with (|209| ) if p 7^ 2 [respectively with 5 > — 1 if p = 2] . 

Remark. Under the hypothesies made on (M, 7) and k, the solution u £ W^^2 s ^^ ^^ equation 

A^u~ku^v (211) 

with V e W^g,_^_2 for some S' such that (5 < (5' < -- + n - 2 is in W^^^ g, ii p > |. Indeed u e Wf+2 5 and k £ T^f ^-^j 
imply that ku £ W^^g,,^^' si'^ce s < 2s + 2 - ^ if p > |, 5" < (5 + (5 + ^). Since u satisfies 

A.,u = ku + v e Vl^,^ i„t(,„,,,)+2 . (212) 

we have 

^^W!^2.JnfiS",S')- (213) 

An induction argument shows that u £ T^f , 2 g/ ■ 
2. Non-Linear Equations 

We suppose that the function / is smooth in y and Wja_|_2 in x. To make it more transparent we take / as a finite 
sum of products of functions of x by functions of y, as it appears in the hamiltonian constraint: 

Q 
f{x,y)=Y,ap{x)bp{y). (214) 

p=o 

We make the following hypothesis: 
Hypothesis H {Wig). 

1. There exists an interval / C M such that the 6's are smooth functions oi y £ I. 

2. The a's are functions on M in M^f ^ , 2- 

Lemma 1. Under the hypothesis H{Wlg) the function on M given by a; i-^ f{x,ip{x)), denoted in the sequel 
/(x, if), has the following properties when ip is continuous and takes its values in a closed interval [i, m] C /: 

I. f{x,^)£WPg^, iis>0. 
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2. If s > 2 , 5 > _H and D(^ e W^I^-i^S'+i ^ith S' > -^ , s > s' > ^ then /(x, <^) G T4^f,^5+2 • 

Proof. Part 1 is trivial. To prove part 2 one uses the calculus derivation formulas and the multiplication properties 
of weighted Sobolev spaces. 

Definitions. A C^ function tp- on M is called a subsolution of A^ip — f{x, ip) if it is such that on M, 

A^(^_>/(a;, </._). (215) 

A C^ function if+ is called a supersolution if on M 

A^ip+<f{x,ip+). (216) 

Theorem 1 (existence). Let (M, 7) be a (p, a, p) asymptotically euclidean manifold cr > - + 2 and f{x,y) a function 
satisfying the hypothesis {H) with s > — , 5 > — - . Suppose the equation A^Lp — f{x, ip) admits a subsolution </?_ 
and a supersolution (p_|_ such that on M: 

i<V-<'P+<m, [£,m]cl (217) 

and 

limip^ < 1 , liiinp+>l . (218) 

00 00 

Suppose that Dip- , Dip-\- G ^f/_i 5/ 1 1 , s' > s ^ 6' > — — . Then the equation admits a solution (f such that 

Tl Ti 

P- <'P<if+ , I- ip e W^^2 5 ^itt^ <5</3and - - < S < n - 2 - - . (219) 



Proof. We construct a solution by induction, starting from ip- 

s,S+2 



Let fc be a positive function on M such that k e W^ g,2 ''^^'^ ^t each point x £ M 



k{x)> sup f'y{x,y). (220) 

£<y<m 

Such a function exists by the hypothesis made on /. 
We set Lp\ = 1 + Ml. The linear elliptic equation for ui: 

A^ui~kui^ f{x,ip^) -k{(p^ -I) (221) 

has one solution ui G T^f_i_2 a ^ C^, since the right hand side is in W^g,2 ^^i^ the operator on the left is injcctive 
from VFf I 2 s ^^^^ ^s 5+2 under the hypothesis made on s and S. The function (pi tends to 1 at infinity. 

We deduce from the equality and the inequalities satisfied respectively by (pi and ip^ the following inequality: 

A^((^i - p_) - k{pi - ^_) < ; (222) 

hence, by the maximum principle lemma, since ipi — p^ tends to c > at infinity, 

pi>p- on M . (223) 

Also, 

A^((/?+ - pi) - k{p+ - pi) < f{x, p+) - f{x, p^) - k{p+ - p^) , (224) 

and 

/(x, p+) - fix, p^) - {p+ - p^) J f'y{x, p^ + t{p+ - p^))dt . (225) 

By the hypothesis made on fc, (/3_|_, and 931 we have on M 

A^{p+ ^ pi) - k{p+ - pi) < , (226) 

hence 
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tpi < </'+ ■ (227) 

The induction formula is, with </?„ = 1 + u„: 

AgUn - kun = fix, ipn-i) - ku^-i ■ (228) 

We suppose that ^p exists for Q < p < n — \ with Lp^ = Lp_ and Up G W^f-i_2 s ^'^'^ 1 < p < n — 1 and that for these p's 

ip- < ipp-i < ipp < ip+ ■ (229) 

The elhptic theory shows u„ S ^^1+2 a exists. The functions tpn are continuous, even C^ since s > -, and tend to 1 
at infinity. The equahty resulting from the equations satisfied by ipp when p < n — 1 gives 

AgCpn-l - kipn-l = fix, (Pn-2) - kipn-2 , (230) 

Agipp-k(pp = fix,Lpp_i)- kipp^i . (231) 

One deduces then from the maximum principle lemma that on M 

(/Sn-l < Vn ■ (232) 

Analogously one uses the maximum principle and the inequality deduced from the equation and inequality satisfied 

by ipn and ip+, 

A^((/3„ - ip+) - kiipn - V+) > fix, </?„-i) - fix, tf+) - kiipn-1 - 'P+) , (233) 

to show that ipn-i < y^+ implies ipn < ^+ ■ 

The sequence of continuous functions (/?„ has been shown to be pointwise increasing and bounded. It is therefore 
converging at each point a; G M to a limit ifix) = 1 + uix), with (p_ < (p < (p+ . 

To show that (p is a solution of the given equation and (p — 1 G W^_^2 s ^*^ proceed as follows. Since the (p„ are 
continuous and take their values in the interval [^, m], the functions f{x, ip„) — kun belong to Wq ^,2 with uniformly 



bounded norms. The linear elliptic inequality (following from (228)) 



shows that the sequence u„ is uniformly bounded in the VF|'^ norm. Since W2 g is compactly embedded in Wf g, for 
any S' < S, there is a subsequence, still denoted u„, which converges in Wf g, norm to a function u g W2 g ■ 

The functions /(a;, (/?„) converge to fix, if) in the W^g,,2 norm because of the inequality, which is satisfied if 
s > li ^' > _» 

\\fix,ip)~ fix,ipn)\\w^ < C||^-(p„||wP IJFilltyP , (235) 

where C depends only on (Af , e) and Fi G ^f 5 i 2 ^^ ^ fnnction on Af , which exists by the hypothesis on /, snch that 

Fiix)> sup |/;(a:,y)|. (236) 

ye[i,m] 

These convergences imply that the limit cp ~ l + u satisfies the equation in a generalized sense. From the linear theory, 
we find that the equation satisfied by u and the fact that (cf. above lemma) that fix,ip) £ Wf g^2 (also € W2 g_^_2 
since u € W2 g) that u G W^ g . An induction argument completes the proof that u G Wf^2 5 ■ 

The theorem holds with s = if / is an increasing function of y. An example is treated in Section [X[ 
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